Transverse magnetization and torque in asymmetrical mesoscopic superconductors 
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We show that asymmetrical mesoscopic superconductors bring new insight into vortex physics 
where we found the remarkable coexistence of long and short vortices. We study an asymmetrical 
mesoscopic sphere, that lacks one of its quadrants, and obtain its three-dimensional vortex patterns 
by solving the Ginzburg-Landau theory. We find that the vortex patterns are asymmetric whose 
effects are clearly visible and detectable in the transverse magnetization and torque. 
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The small volume to surface ratio in mesoscopic su- 
perconductors makes their properties distinct from the 
bulk materials.^ For instance, in the absence of pinning 
and anisotropy, mesoscopic superconductors have an in- 
trinsic magnetic hysteresis whereas bulk samples do not. 
Mesoscopic superconductors are intrinsically metastable 
since they can be in exited vortex states due to an ener- 
getic barrier upholding the decay to the groundstate^^A^ 
Several experimental techniques have been developed to 
detect mesoscopic vortex states^^^^^ and to reveal their 
new and interesting physical properties ^^i^ Here we ob- 
tain the transverse magnetization and torque that result 
from the asymmetrical coexistence of vortices with differ- 
ent lengths in mesoscopic superconductors. The general 
aspects of our results apply to irregularly shaped grains 
and, consequently, can also be relevant to the under- 
standing of inhomogeneous bulk samples with granular 
struct urei^ii^iii 

In this paper we obtain and compare the vortex pat- 
terns of two similar but distinct mesoscopic systems, 
namely, a full sphere with radius R = 6.0^ and an 
asymmetrical sphere, with the same radius but lacking 
one of its quadrants. Hereafter we call it minus- quarter 
sphere and present here its transverse magnetization and 
torque. There are many reasons for the magnetic mo- 
ment of a sample with volume = VM^ not to be 
aligned along the applied magnetic field, H. A direct 
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of the torque 
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brings insight about these distinct contributions ^li^ 
where is the angle between H and M. For example, 
vortices pinned by point-like normal inclusions (defects) 
are major contributors to torque due to their inability 
to follow a rotating applied field. This contribution is 
usually avoided as one seeks out true information about 
intrinsic properties of the superconducting state. The 
sample's shape and size is also another contributor and 
provides the basis for very high precision magnetization 
measurements using torque. Basically it relies on the 
thin film limit, where vortices remain straight and single 
oriented inside the sample upon field rotation, resulting 
in a magnetic sign always perpendicular to the major 



(g) 



(h) 



(d) 



(e) 



(a) 



(b) 




FIG. 1: (Color online) Three-dimensional iso-density plots 
taken at 20% of the maximum density Each iso-contour 
is a single surface made of the vortices and thesuperconduc- 
tor surface. The magnetic field is directed from bottom-to-top 
with H/Hc2 = 0.16 (a), 0.29 (b), 0.35 (c), 0.27 (d), 0.39 (e), 
0.49 (f), 0.35 (g), 0.48 (h), and 0.54 (i). The figures cor- 
respond to the full sphere (a)-(c) and to the minus- quarter 
sphere, positioned in two ways, ((d)-(f) and (g)-(i), with re- 
spect to the field. 



surfaces. The general misalignment between the applied 
field and the sample's major axis produces a torque, even 
for a regularly shaped superconductor. Hence to obtain a 
torque it suffices that the circulating supercurrent yields 
a magnetic field response not oriented along the applied 
field that gave rise to it. Asymmetry leads to such effect 
and we show here that it provides a useful tool to probe 
the physics of mesoscopic superconductors. 

In this paper we provide the first study of a meso- 
scopic superconductor with a transverse (Mx) magneti- 
zation, thanks to the use of a truly three-dimensional 
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approach. ^^'^^ The found effects in the transverse magne- 
tization and in the torque are within reach of present ex- 
perimental detection. For instance, a coherence length 
of ^ :^ 100 nm leads to Hc2 = ^o/^tt^^ 300 Gauss. 
If the penetration length is comparable to the sphere's 
radius R = 0.6 /im, then the transverse magnetization is 
approximately ten times smaller, — 10~^Hc2 — 0, 3 
Gauss than the longitudinal one Mz — 10~^Hc2 — 3 
Gauss. The magnetic moment of this sphere is /i ~ 
MzV - 7.0 X 10-^2 emu (V = 727t^^ 7.63 x 10"^^ 
cm^ for the minus- quarter sphere). The torque, Ty = 
-M^HzV lO-^H'^^y becomes Ty = 4.0 x 10-^^ erg. 

Fig. [T] summarizes the novel aspects brought by asym- 
metry showing the Cooper pair density, |?/^(r)p, obtained 
from our numerical simulations for selected values and 
orientations of the applied field. States with at most 
three vortices were selected to help to understand the 
basic effects caused by asymmetry. For all cases the 
magnetic field is oriented along the bottom-to-top direc- 
tion. In those mesoscopic systems vortices are naturally 
curvectf^ because the entrance and exit of the vortex lines 
must be perpendicular to the surface, ^^«« a fact that causes 
a closely packed configuration near to the equatorial re- 
gion, as previously found. 

Fig. [IJd) best represents the effects of the asymme- 
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FIG. 2: (Color online) The magnetization for a full field cycle 
is shown here, (a) for the full and (b)-(c) for the minus- 
quarter sphere. The ascending (red) and descending (blue) 
branches are labeled for the first nine vortex states. Selected 
iso-plots obtained in ascending and descending fields show the 
degeneracy of the 5 state. The full non-intersecting magneti- 
zation branches are shown just for the full sphere. 
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FIG. 3: (Color online) The ascending (red) and descending 
(blue) torque branches are shown here in (a) sub-plot and the 
corresponding misalignment angle in (b) sub-plot. 



try where theposition of the short vortex very close to 
the center of the minus- quarter sphere is determined 
by: (i) the spherical shell pushing the vortex near the 
central axis where the shielding current offers maximal 
stability^^ and, (ii) the fiat surface of the missing quad- 
rant where the vortex length is half of the one for a full 
sphere. As a consequence, the vortex is displaced from 
the central axis, breaking the rotational symmetry in the 
northern hemisphere resulting in transverse magnetiza- 
tion and torque. The second vortex has the same faith 
of the first one, it is also a short vortex that gets as close 
as possible to the central axis (Fig. [TJe)). The situation 
changes for the third vortex, possibly because R = 6.0^ 
is just too small to accommodate it in the quadrant, forc- 
ing its nucleation in the full hemisphere where it becomes 
the first long vortex. Asymmetry brings peculiar, not al- 
ways intuitive, features. For instance the rotation of the 
minus- quarter sphere by 45° leads to new vortex config- 
urations not accessible in case of the previous alignment. 
Now the long vortices nucleate first and short vortices 
enter in pairs afterwards, as shown in Figs.[T](g)-(i), with 
increasing field. In summary asymmetry brings new sit- 
uations that can bring further information about vor- 
tex patterns in mesoscopic superconductor. They are de- 
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tectable in the transverse magnetization and the torque, 
bringing inside information about vortex patterns in such 
systems. 

The Ginzburg-Landau theory, derived near the critical 
temperature, is apphcable at much lower temperatures 
for the mesoscopic superconductors, where it was found 
to give a fair description of its magnetic properties i^ii^ 
We solve the full three-dimensional Ginzbug-Landau the- 
ory for the asymmetrical superconductor by the numer- 
ical procedure described in Refs. 17.23, that minimizes 
the free energy with the boundary conditions included. 
In this approach the superconductor shape is character- 
ized by a step-like function in the free energy (1 inside 
the superconductor and outside). Here we also in- 
clude an intermediate step inside the superconductor (0.8 
within a layer of thickness 0.5^) to stabilize the states for 
fields near i^c3, according to Ref. [l^. The present results 
are under the approximation that the applied field uni- 
formly permeates the mesoscopic superconductor so that 
the magnetic shielding is safely ignored. A discretized 
gauge invariant version of the Ginzburg-Landau free en- 
ergy is minimized in a cubic mesh with 46 x 46 x 46 
points^^ and the magnetization is directly determined 
from M = const (l/2c) J r x J dv/V^ where J is the 
supercurrent in the mesoscopic superconductor volume 
and const a free parameter associated to the demagne- 
tization tensor^^^^^ and so, obtained in our numer- 
ical simulations from the zero vortex branch (Meissner 
phase) at small H: H ^ AttD • M = 0. We apply the in- 
verted relations Mx = {—DzzHx-\-DxzHz)/det and = 
{DzxHx — DxxHz) / det with det = DzzDxx — DxzDzx in 
our numerical analysis to find that for the full sphere 
the demagnetization tensor reduces to a single scalar, 
as expected, which is adjusted to the well known 
value^iiS^ of I) = 1/3 and for the minus- quarter sphere 
Dxx — Dzz^ Dxz = Dzx^ and also that Dzz/D = 1.20 
and Dzx/D = —0.184. In the rest of the paper we use 
the gaussian system of units and express both the mag- 
netization and the applied field in terms of i^c2, the bulk 
upper critical field, and the torque per volume in units 

Fig. [2] shows the magnetization of the full and 
the minus- quarter spheres. We use the notation 
long -\- short {long^ short) to express the minus-quarter 
sphere states, whose first nine ones, labeled in Figs. [2jb) 
and (c), are given by: (0,0), 1 (0,1), 2 (0,2), 
3 (1,2), 4 ^ (1,3), 5d (1,4), 5a ^ (2,3), 6 ^ 
(2, 4), and 7 (2, 5). Notice that there are two 5 vortex 
states, 5a, and 5(i, the first being the lowest in energy, 
and therefore reached by the ascending field. These two 
states can coexist under the same field although Fig. [2] 
does not show the corresponding prolongation of the 5a 
and 5d curves. This degeneracy is a new and interest- 
ing aspect brought by asymmetry which shows that the 
total number of vortices not always determines iniquely 
the magnetization branch. We found 24 states for the full 
spheres and at least 26 for the minus- quarter sphere. The 
higher vortex states, which live near to the normal state. 



become increasingly difficult to identify and so, required 
more precision than considered here. The magnetization 
curves were obtained to a full cycle field sweep that starts 
at zero, reaches the maximum possible value (Hcs) (as- 
cending branch, red color), and finally is lowered back 
to zero (descending branch, blue color). Figs, and 
(b) show the magnetization component along the field, 
M^, with the expected feature of rounded branches for 
ascending field and of a tooth-saw structure for descend- 
ing field. This behavior has been extensively discussed 
both theoretically^^i2ii2^ and also experimentally.^ The 
full stability region for the different vortex configuration 
is only displayed here for the full sphere (Fig. [2ja)). Its 
construction requires numerous sweeps starting at inter- 
mediate applied field values. Notice the distinct features 
of Mx in comparison to the M^. The ascending and de- 
scending branches are entangled and have an underlying 
structure that approximately breaks in distinct families, 
associated to odd and even number of vortices. These 
branches display a wiggled structure due to vortex repo- 
sitioning that cause abrupt changes in vortex length, and 
so in energy and in current fiow. 

The vortex entrance or exit can be detected through 
the transverse magnetization, whose sign changes accord- 
ingly, as seen in Fig. [21(c). Mesoscopic superconductors 
are known to be able to exhibit a paramagnetic sig- 
nal when in a metastable vortex state^^'^^ but this is 
not the case for the sphere, as the dominant component 
{Mz) remains steadily oppositely oriented to the applied 
field. However a change of sign in Mx implies that the 
torque fiips direction, as shown in Fig. [3l Moreover the 
torque is proportional to the applied field and it gets 
naturally enhanced in the regime between Hc2 and i^cs, 
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FIG. 4: (Color online) The torque is shown upon rotation of 
the fixed applied field {H — 0.40) in 5° increments. The inset 
shows the angle between M and H versus the applied field 
angle. Selected iso-plot configurations for 110°, 130°, 150°, 
and 180° are also displayed here. 
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contrary to the magnetization, which gets dim in this re- 
gion. The magnetization angle, 180° — ^, shows that the 
misahgnment with the apphed field is about 10° for small 
fields and it oscillates above i^c2- Thus, torque measure- 
ment can be very helpful in the detection of high vortex 
states where the onset of a vortex configuration causes 
an abrupt change of sign. As an example. Fig. [H shows 
the torque for a rotating field of intensity H = 0.4. The 
selected initial configuration for 0° is the ground state for 
this field, made of two short vortices (see Fig. [Tl|e)). As 
the field rotates by an angular increments of 5°, the two 
short vortices follow along and this configuration remains 
stable until 110°, albeit the fact that during this process 
the torque fiips sign four times! The torque suffers an 
abrupt drop at 115° signaling the entrance of a new vor- 
tex pattern. One of the short vortices is replaced by a 
long vortex and the remaining short vortex moves to the 
center of the minus- quarter sphere. This long and short 
vortex configuration, i.e. (1,1) lasts until the next torque 
transition at 130°, where a single long vortex becomes the 



accessed configuration (from 135° until 150°). Thus the 
Fig. [11(g) state is reached in this way. The last abrupt 
drop brings back the (1, 1) configuration from 155° until 
180°. The final state does not coincide with the initial 
one, and is not included as one of the available stable 
configurations listed in Fig. [TJ Thus it is a very high ex- 
cited state much above the ground state. Nevertheless it 
became accessible by virtue of the torque. 

In conclusion we have shown here that the extra com- 
plexity of vortex patterns in asymmetrical mesoscopic 
superconductors, under a tilted field, can be used as a 
tool for vortex detection through the mesurement of the 
transverse magnetization and the torque. 

A. R. de C. Romaguera acknowledges the Brazilian 
agency CNPq for financial support. M. M. Doria thanks 
CNPq, FAPERJ, Instituto do Milenio de Nanotecnolo- 
gia (Brazil) and BOF/UA (Belgium). F. M. Peeters ac- 
knowledges support from the Flemish Science Foundation 
(FWO-Vl), the Belgian Science Policy (lUAP) and the 
ESF-AQDJJ network. 



* Electronic address: ton@if.ufrj.br 

^ A. T. Dorsey, Nature (London) 408, 783 (2000). 

^ B. J. Baelus, F. M. Peeters, and V. A. Schweigert, Phys. 

Rev. B 63, 144517 (2001). 
^ P. S. Deo, V. A. Schweigert, and F. M. Peeters, Phys. Rev. 

B 59, 6039 (1999). 
^ A. K. Geim, I. V. Grigorieva, S. V. Dubonos, J. G. S. Lok, 

J. C. Maan, A. E. Filippov, and F. M. Peeters, Nature 

(London) 390, 259 (1997). 
^ C. A. Bolle, V. Aksyuk, F. Pardo, P. L. Gammel, E. Zeldov, 

E. Bucher, R. Boie, D. J. Bishop, and D. R. Nelson, Nature 

(London) 399, 43 (1999). 
^ A. K. Geim, S. V. Dubonos, J. J. Palacios, L V. Grigorieva, 

M. Henini, and J. J. Schermer, Phys. Rev. Lett. 85, 1528 

(2000). 

^ L V. Grigorieva, W. Escoffier, J. Richardson, L. Y. Vin- 
nikov, S. Dubonos, and V. Oboznov, Phys. Rev. Lett. 96, 
077005 (2006). 

^ A. Kanda, B. J. Baelus, F. M. Peeters, K. Kadowaki, and 

Y. Ootuka, Phys. Rev. Lett. 93, 257002 (2004). 
^ A. R. Jurelo, J. V. Kunzler, J. Schaf, P. Pureur, and 

J. Rosenblatt, Phys. Rev. B 56, 14815 (1997). 

C. Howald, P. Fournier, and A. Kapitulnik, Phys. Rev. B 

64, 100504 (2001). 

T. Ekinoa et al, Physica C 426-431, 230 (2005). 
^2 M. M. Doria and L G. de Oliveira, Phys. Rev. B 49, 6205 
(1994). 

V. G. Kogan, Phys. Rev. Lett. 89, 237005 (2002). 

Y. Wang, L. Li, M. J. Naughton, G. D. Gu, S. Uchida, and 

N. P. Ong, Phys. Rev. Lett. 95, 247002 (2005). 



D. E. Farrell, C. M. Williams, S. A. Wolf, N. P. Bansal, 
and V. G. Kogan, Phys. Rev. Lett. 61, 2805 (1988). 

H. Won and K. Maki, Europhys. Lett. 34, 453 (1996). 
M. M. Doria, A. R. de C. Romaguera, and F. M. Peeters, 
Phys. Rev. B 75, 064505 (2007). 

A. R. de C. Romaguera, M. M. Doria, and F. M. Peeters, 
cond-mat/0609576 (2006). 

P. S. Deo, V. A. Schweigert, F. M. Peeters, and A. K. 
Geim, Phys. Rev. Lett. 79, 4653 (1997). 

A. R. de C. Romaguera, M. M. Doria, and F. M. Peeters, 
Phys. Rev. B In press (2007). 

E. H. Brandt, J. Low Temp. Phys. 42, 557 (1981). 

B. J. Baelus, D. Sun, and F. M. Peeters, Phys. Rev. B In 
press (2007). 

M. M. Doria and A. R. de C. Romaguera, Europhys. Lett. 
67, 446 (2004). 

J. A. Osborn, Phys. Rev. 67, 351 (1945). 

M. Beleggia, M. D. Graef, and Y. T. Millev, J. Phys. D: 

Appl. Phys. 39, 891 (2006). 

V. A. Schweigert, F. M. Peeters, and P. S. Deo, Phys. Rev. 
Lett. 81, 2783 (1998). 
^'^ V. A. Schweigert and F. M. Peeters, Phys. Rev. Lett. 83, 
2409 (1999). 

B. J. Baelus and F. M. Peeters, Phys. Rev. B 65, 104515 
(2002). 

J. J. Palacios, Phys. Rev. Lett. 84, 1796 (2000). 

A. K. Geim, S. V. Dubonos, J. G. S. Lok, M. Henini, and 

J. C. Maan, Nature (London) 396, 144 (1998). 



